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Abstract 

Conventional thermo-statistics address infinite homogeneous systems within the 
canonical ensemble. However, some 170 years ago the original motivation of ther- 
modynamics was the description of steam engines, i.e. boiling water. Its essential 
physics is the separation of the gas phase from the liquid. Of course, boiling water 
is inhomogeneous and as such cannot be treated by conventional thermo-statistics. 
Then it is not astonishing, that a phase transition of first order is signaled canoni- 
cally by a Yang-Lee singularity. Thus it is only treated correctly by microcanonical 
Boltzmann-Planck statistics. This was elaborated in the talk presented at this con- 
ference. It turns out that the Boltzmann-Planck statistics is much richer and gives 
fundamental insight into statistical mechanics and especially into entropy. This can 
be done to a far extend rigorously and analytically. The deep and essential differ- 
ence between "extensive" and "intensive" control parameters, i.e. microcanonical 
and canonical statistics, was exemplified by rotating, self-gravitating systems. In 
the present paper the necessary appearance of a convex entropy S{E) and the 
negative heat capacity at phase separation in small as well macroscopic systems 
independently of the range of the force is pointed out. 
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1 Introduction 



Since the beginning of thermodynamics in the first half of the 19. century its 
original motivation was the description of steam engines and the liquid to 
gas transition of water. Here water becomes inhomogeneous and develops a 
separation of the gas phase from the liquid, i.e. water boils. 

A little later statistical mechanics was developed by Bottzmann(0) to explain 
the microscopic mechanical basis of thermodynamics. Up to now it is gener- 
ally believed that this is given by the Boltzmann-Gibbs canonical statistics. 
As traditional canonical statistics works only for homogeneous, infinite sys- 
tems, phase separations remain outside of standard Boltzmann-Gibbs thermo- 
statistics, which, consequently, signal phase-transitions of first order by Yang- 
Lee singularities. 



It is amusing that this fact that is essential for the original purpose of thermo- 
dynamics to describe steam engines was never treated completely in the past 
150 years. The system must be somewhat artificially split into (still macro- 
scopic and homogeneous) pieces of each individual phase (0). The most in- 
teresting configurations of two coexisting phases cannot be described by a 
single canonical ensemble. Important inter-phase fluctuations remain outside, 
etc. This is all hidden due to the restriction to homogeneous systems in the 
thermodynamic limit. 



Also the second law can rigorously be formulated only microcanonically: Al- 
ready Clausius (J3|) distinguished between external and internal entropy gener- 
ating mechanisms. The second law is only related to the latter mechanism (jij), 
the internal entropy generation. Again, canonical Boltzmann-Gibbs statistics 
is insensitive to this important difference. 



For this purpose, and also to describe small systems like fragmenting nuclei 
or non-extensive ones like macroscopic systems at phase-separation, or even 
very large, self-gravitating, systems, we need a new and deeper definition of 
statistical mechanics and as the heart of it: of entropy. For this purpose it is 
crucial to avoid the thermodynamic limit. 



As the main aspects of this new thermodynamics were published in 
0: S 0) I skip here to repeat all the arguments. Instead I will stress here 
only the fact that negative heat capacity and convex entropy can be seen at 
proper phase transitions of 1. order, i.e. at phase separation, in small as well 
in macroscopic systems independently whether they have long or short range 
interactions. As there was a hot discussion at this conference about this point, 
it seems necessary to repeat the arguments here. 
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2 Negative heat capacity at phase-separation can also be seen in 
macroscopic systems independently of the range of the interac- 
tion. 



The argument is simple c.f.(9): At phase separation the weig ht e s ^- E ' T of 
the configurations with energy E in the canonical partition sum 

oo 

Z(T) = J e s ^- E ' T dE (1) 



becomes bimodal, at the transition temperature it has two peaks, the "liquid" 
and the "gas" configurations which are separated in energy by the latent 
heat. Consequently S(E) must be convex (d 2 S/dE 2 > 0, like y = x 2 ) and 
the weight in (1) has a minimum at E m i n between the two pure phases. Of 
course, the minimum can only be seen in the microcanonical ensemble where 
the energy is controlled and its fluctuations forbidden. Otherwise, the system 
would fluctuate between the two pure phases by an, for macroscopic systems 
even macroscopic, energy AE ~ E\ at oc iV of the order of the latent heat 
in clear contrast to the usual assumption of the fluctuations in the canonical 
ensemble 5E oc \fN . The heat capacity is 



C v (E min ) = dE/dT = - 



(dS/dEf 



d 2 S/dE 2 



< 0. 



(2) 



I.e. the convexity of S(E) and the negative microcanonical heat capacity are 
the generic and necessary signals of any phase-separatio ml u\) . 

This "convex intruder" in S(E) with the depth AS sur f(E min ) has a direct 
physical significance: Its depth is the surface entropy due to constraints by the 
existence of the inter-phase boundary between the droplets of the condensed 
phase and the gas phase and the corresponding correlation. AS sur f(E min ) is 
directly related to the surface tension per surface atom (with number N sur f) 
of the droplets. 

0~surfl J-tr — ^7 l<Jj 



N, 



surf 



In my paper together with M.Madjet (jllh we have compared the values of 
AS sur f(E min ) calculated by Monte-Carlo using a realistic short range interac- 
tion with the values of the surface tension of the corresponding macroscopic 
system. In these calculations we used the empirical liquid drop parameters for 
the ground-states energies of the different clusters as given by dJJ). Table (1) 
shows the scaling behavior of AS sur f(E min ) with the size iV of the system. 
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Fig. 1. Microcanonical Monte-Carlo (MMMC) (fill : fioh simulation of the entropy 
s(e) per atom (e in eV per atom) of a system of N = 1000 sodium atoms at an ex- 
ternal pressure of 1 atm. At the energy e < e\ the system is in the pure liquid phase 
and at e > e3 in the pure gas phase, of course with fluctuations. The latent heat per 
atom is qi at = e% — e\. Attention: the curve s(e) is artificially sheared by subtracting 
a linear function 25 + e * 11.5 in order to make the convex intruder visible. s(e) is 
always a steep monotonic rising function. We clearly see the global concave (down- 
wards bending) nature of s(e) and its convex intruder. Its depth is the entropy loss 
due to additional correlations by the interfaces. It scales oc N^ 1 ^ 3 . From this one 
can calculate the surface tension per surface atom cr sur f/Tt r = As sur f * N/N sur f. 
The double tangent (Gibbs construction) is the concave hull of s(e). Its derivative 
gives the Maxwell line in the caloric curve e(T) at Tt r . In the thermodynamic limit 
the intruder in s(e) would disappear and s(e) would approach the double tangent 
from below, not of course in S(E), which remains deeply convex: The probability of 
configurations with phase-separations is suppressed by the (infinitesimal small) fac- 
tor e~ N2 3 relative to the pure phases and the distribution remains strictly bimodal 
in the canonical ensemble in which the region e\ < e < e^ of phase separation gets 
lost. 

Roughly AS surf (E) oc N 2 ' 3 and one may argue that this will vanish compared 
to the ordinary leading volume term S vo i(E) oc N. However, this is not so as 
S vo i(E) at energies inside the phase-separation region (the convex intruder) 
is the concave hull of S(E) (its slope gives the Maxwell construction of the 
caloric curve T(E)). It is a straight line and its curvature d 2 S vo i/dE 2 = 0. 
Consequently for large N 

d 2 S/dE 2 ~ d 2 S vol /dE 2 + d 2 AS surf /dE 2 + • ■ ■ x d 2 AS surf /dE 2 (4) 

and the depth of the intruder AS sur f(E min ) = a/T tr * N sur j ~ jV 2 / 3 goes 
to infinity in the thermodynamic limit. Of course, the ubiquitous phenomena 
of phase separation exist only by this reason. It determines the (negative) 
heat capacity as in eq.(2). The physical (quite surprising) consequences are 
discussed in (J^l 0) . 
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Table 1 

Parameters of the liquid-gas transition of small sodium clusters (MMMC- 
calculation (fill: Hob) in comparison with the bulk for a rising number N of atoms, 

2/3 

N sur f is the average number of surface atoms (estimated here as Yl ^cluster) °^ an 
clusters with iVj > 2 together. a/T tr = As sur f * N/N sur f corresponds to the surface 
tension. Its bulk value is adjusted to agree with the experimental values of the a s 
parameter which we used in the liquid-drop formula for the binding energies of small 
clusters, c.f. Brechignac et al. (|l2h . and which are used in this calculation |l(ih for 
the individual clusters. 





N 


200 


1000 


3000 


bulk 




Ttr [K] 


940 


990 


1095 


1156 




qiat [eV] 


0.82 


0.91 


0.94 


0.923 


Na 


Sboil 


10.1 


10.7 


9.9 


9.267 






0.55 


0.56 


0.44 






N t 

1 v surj 


39.94 


98.53 


186.6 


oo 




cr/Tt r 


2.75 


5.68 


7.07 


7.41 



Discussions with St.Ruffo are gratefully acknowledged. 
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